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One of the major purposes in studying heavy-ion collisions is to understand the physics of particle interactions at 
both high particle and high energy densities and its dependence on the colliding energy. Thus, based on the free 
£f~) • space cross section (or mean- free-path) arguments, it was suggested long time ago [1] that the freeze-out sequence 
is: photons, K + 's, protons and pions. The transport model calculations then confirmed this expectation except for 
protons which appear to be emitted latest in the mean. This is demonstrated in Fig. 1 where we present the particle 
yields as functions of the emission time in the reaction c.m.s. for Pb + Pb central collisions at 158 A-GeV simulated 
$_i ' with the RQMD (v2.3) code. Nowadays, in connection with the experiments at SPS with heaviest nuclei and the 
future experiments at RHIC and LHC, there are of special interest the effects of the eventual phase transitions, like 
the production of quark-gluon plasma or strangeness distillation, leading to specific delays in particle production. 
Therefore the measurement of the freeze-out sequence may be able to provide crucial information on the properties 
of the fireball produced in heavy-ion collisions. 

In this report, we address the problem of relative delays and spatial shifts among unlike particles. We use the 
i-H transport RQMD (v2.3) code to generate the 158 A-GeV Pb+Pb central collisions. The correlation functions for 

unlike particle pairs: (tt + ,K + ), (tt + ,p) and (K + ,p) are then calculated with the final state interaction code [2,3]. 
(N : In order to minimize the effect of a fast longitudinal motion of the particle sources, we consider here only the mid- 
rapidity region \y\ < 1 and calculate the two-particle correlation functions in the longitudinally co- moving system 
(LCMS). In this system the particle pairs are emitted transverse (x = out) to the reaction axis (z — long). In Fig. 2 
we present the distributions of the relative ^-coordinates of the emission points and of the relative emission times 
, calculated in RQMD. The distributions of the relative y- and z-coordinates show practically no asymmetry as expected 
from symmetry reasons (the azimuthal symmetry and the symmetry of the initial system together with the choice of 
"*^h | a symmetric mid-rapidity window). 

As pointed out in [4] , we can access the information about these asymmetries exploiting the fact that the final state 
interaction, determining the correlation function of two non-identical particles, depends on the orientation of the 
relative 4-coordinates of the emission points through the scalar product q-r* , where r* — fj — f% = {Ax*, Ay*, Az*} 
and q = pi = —p^ are the relative 3-coordinates of the emission points and half the relative momentum in the 
two-particle c.m.s. Since this system moves in the x-direction with the velocity v± = \pi±_ +Pi±\/(Ei + E2) and the 
Lorentz factor yj_ = (1 — v 2 ^)^ 1 / 2 , we have: 
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Ax* = 7_l(Ax - v±At), Ay* = Ay, Az* = Az. (1) 

Since further (q- r*) = q x (Ax*) + q y (Ay*) +q z (Az*) and, in our case Ay = Az = 0, we can search for the LCMS space 
(out) and time asymmetries with the help of the correlation functions TZ+ and 1Z- corresponding to q x = q-Vj_/v± > 
and < 0, respectively [4]. In the case when the time differences dominate over the spatial ones, i.e. when v\t\ r 
on average, the ratio of the correlation functions 1Z+ and 7£_ directly measures the mean relative emission time, 
including its sign [4] . Noting that for particles of equal masses the sign of the scalar product q ■ v coincides with the 
sign of the velocity difference v% — V2 , we can see the simple classical meaning of the above selection. It corresponds 
to the intuitive expectation of different particle interaction in the case when the faster particle is emitted earlier as 
compared to the case of its later emission. In former case the interaction between the two particles will be weaker 
and the correlation function 1Z + will be closer to unity than TZ- . 

To clarify the origin of the asymmetry effect, we will follow the arguments given in [4]. On the usual assumptions of 
a small particle phase-space density and sufficiently smooth behavior of the single-particle spectra, the two-particle 
correlation function is determined by the modulus squared of the two-particle amplitude ip_g(f**) averaged over the 
relative coordinates of the emission points in the two-particle c.m.s. [2]: 

R(Pun) = <IV>-#*)I 2 >- (2) 

For charged particles the final state interaction is often dominated by the Coulomb interaction so that 
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where £ = gr* + gr* = qr* (1 + cos #*), 77 = l/(k*a), a = l/(^ziZ2e 2 ) is the Bohr radius (p is the reduced mass of the 
two particles, z\z 2 e 2 is the product of the particle electric charges; a = 248.5, 222.5 and 83.6 fm for the ir + p, K + p 
and 7r + K + systems, respectively), <5g = argr(l + it]) is the Coulomb s-wave phase shift, 

A c (r}) = 27r?7[cxp(27n7) - l]" 1 (4) 

is the Coulomb factor - modulus squared of the non-relativistic Coulomb wave function at zero distance - introduced 
already by Fermi in his theory of /3-decay, 

F(a, l,z) = l + az + a(a + l)z 2 /2! 2 + a(a + l)(a + 2)z 3 /3! 2 + • • • (5) 

is the confluent hypergeometric function. At small relative momenta qr* < 1 and r* <C |a|, we can then write 

R( Pl ,p 2 ) = A c (r,) [1 + 2(r*(l + cos0*))/a + • • •] 

= A c (r]) [l + 2(r*)/a + 2cosip- (Ax*)/a + 2sinip ■ (sin0- (Ay*) + cos0 • (Az*))/a + • • •] , (6) 

where ip and <j> are the polar and azimuthal angles of the vector q with respect to the velocity vector (x-axis): 
q = qjcos ■0, sin -0 sin 0, sin^cos^} and cosi/> = qv±/(qv±). Noting that for particles with the same (opposite) 
charges the Bohr "radius" a > (< 0) and the Coulomb factor A c (rf) approaches unity from below (above), we 
can immediately see that the interaction (and thus the correlation) is stronger if (qr*) = (qr* cos 9*) < 0. In the 
classical limit of qr* » 1 this corresponds to the expectation of stronger interaction between particles moving after 
the emission (on average) towards each other in their c.m.s. 1 

Averaging in Eq. (6) over nearly isotropic distribution of the vector q at q — ► (nearly uniform cosip- and 
(^-distributions) or recalling that in the considered case (Ay*) = (Az*) = 0, we get for the ratio 1Z+/1Z- at q — > 

TZ+/U- « 1 + 2((cosV>) + - (cos^}_) (Ax*) /a, (7) 

where the ± subscripts correspond to the averaging over positive/negative values of q x = q cos ijj. The result of this 
averaging depends on the selected kinematic region. For example, if only pairs with the relative momentum q parallel 
or antiparallcl to the pair velocity are selected, then (cos - 0}± = ±1 and the measured asymmetry would be maximal 
[4]. In the case of averaging over a full range of the uniform cos V'-distribution, we have (cos tp)± = ±1/2 and twice less 
the asymmetry. For our kinematical selection, the corresponding angular factor in brackets in Eq. (7) is somewhat 
higher than unity. 

Regarding the effect of the strong final state interaction, since the magnitude of the two-particle scattering amplitude 
/ is usually much less than the mean r* , it only slightly modifies the correlation functions 1Z+, 1Z- and their 
ratio at small q. In particular, Eq. (7) for the ratio 1Z+/7Z- at q — > is modified by the substitution (Ax*) — > 
(Ax*) + Re/ • (Ax* /r*). This means that for neutral particles (when \a\ — > oo) the ratio 1Z+/TZ- — > 1 at q — > 0. 2 

The correlation functions 1Z+, and their ratio are plotted in Fig. 3. We can see that for n + p and tt + K + systems 
these ratios are less than unity at small values of q, while for K + p system the ratio 1Z + /1Z~ practically coincides with 
unity. These results well agree with the mean values of At, Ax and Ax* presented in Table 1. We may see from Eqs. 
(1) and (7) that the absence of the effect in the 1Z+/1Z- ratio for the K + p system is due to practically the complete 
compensation of the space and time asymmetries leading to Ax* w 0. For ir + p system the effect is determined mainly 
by the x-asymmetry. For -k + K + system both the x- and time-asymmetries contribute in the same direction, the latter 
contribution being somewhat larger. 

We have thus demonstrated that the unlike particle correlations allow to directly access the information about 
space-time asymmetries of particle production in ultrarelativistic heavy ion collisions even at present energies. The 
separation of the relative time delays from the spatial asymmetry is, in principle, possible (see Eq. (1)) by studying 
the ratio 1Z+ /TZ- in different intervals of the pair velocity. However, such a study can hardly be done in a model 
independent way due to the velocity dependence of the out asymmetry Ax (particularly, both Ax and v±At vanish 
at zero LCMS pair velocity v±). 



x It should be noted that the two-particle amplitude tp-^(r*) is often misidentified with the usual two-particle wave function 
in the scattering problem ip$(r*) (see, e.g., [5]). This misidentification would obviously lead to the opposite asymmetry effect 
in the correlation function. 

2 In [4] it was not clearly stated that the asymmetry vanishes at q — > only on the absence of the Coulomb interaction. 
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TABLE I. Mean values of the relative space-time coordinates in LCMS (in fm) calculated from RQMD (v2.3) for 158 A-GeV Pb + Pb 
central collisions. 
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pion 20.17 13.76 

Kaon 18.06 12.33 

proton 26.01 12.83 
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FIG. 1. Emission profiles for mid-rapidity (\y\ < 1) positive pions, kaons and protons from 158 A-GeV/c Pb+Pb central 
collisions. 
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FIG. 2. Space-time distributions for mid-rapidity particle pairs ir + — p, ir + — K + , and K + — p. 
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FIG. 3. Unlike particle correlation functions for mid-rapidity particle pairs ir + — p, 7v + — K + , and K + — p. 
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